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Abstract 

A number of practically relevant problems involving the impulsive motion or the rapid 
rotation of bodies immersed in fluid are susceptible to vortex-like instability modes. 
Depending upon the configuration of any particular problem the stability properties of 
any high-wavenumber vortices can take on one of two distinct forms. One of these is akin 
to the structure of Gortler vortices in boundary layer flows whilst the other is similar to 
the situation for classical Taylor vortices. 

Both the Gortler and Taylor problems have been extensively studied when crossflow 
effects are excluded from the underlying base flows. Recently, studies have been made 
concerning the influence of crossflow on Gortler modes and here we use a linearised 
stability analysis to examine crossflow properties for the Taylor mode. This work allows 
us to identify the most unstable vortex as the crossflow component increases and we show 
how, like the Gortler case, only a very small crossflow component is required in order to 
completely stabilise the flow. Our investigation forms the basis for an extension to the 
nonlinear problem and is of potential applicability to a range of pertinent flows. 
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No. NASl-19480 while the author was in residence at the Institute for Computer Applications in Science and 
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Crossflow effects on Taylor vortices 


1. Introduction 

The problem of predicting vortex breakdown is of great importance within the field 
of transition research. In many practical situations it is desirable to influence breakdown 
characteristics: for example in aerospace design the requirement is that transition be 
delayed as long as possible in order to reduce drag whereas in jet intakes turbulence 
is promoted in order to improve fuel mixing. In numerous applications where flow 
curvature is present the stability problem is closely related to structures describing Taylor 
and/or Gortler modes at large Taylor /Gortler numbers. Undoubtedly in many cases the 
experimental flow is fully three-dimensional in character and until recently little theoretical 
work had been conducted on such problems. 

The linear stability of short wavelength Taylor-Gortler vortices was examined by Hall 
[1] for both fully developed and boundary layer flows. The model problem considered was 
the flow contained between a pair of concentric cylinders whose separation is much smaller 
than the radii of the cylinders. A WKB type analysis shows that in the Dean problem (in 
which both cylinders are at rest and the fluid is driven by a constant azimuthal pressure 
gradient) short wavelength vortices are constrained to lie in a thin layer located at a 
known position between the cylinders. In contrast, for the Taylor problem (in which the 
outer cylinder is at rest and the inner cylinder rotates uniformly) the vortex structure is 
trapped at the inner cylinder wall. These two idealised cases capture the essence of a large 
number of important practical flows. For example, the Dean structure forms the basis for 
describing the evolution of Gortler vortices in spatially varying boundary layers (although 
extra non-parallel effects do need to be considered) whereas the wall-bounded vortices are 
relevant to a number of impulsively started spin-up problems, [2]. The linear stability 
theories for the Taylor and Dean cases have been extended to weakly and fully nonlinear 
accounts, [3-6]. 

The first work concerned with the effect of adding a crossflow component to the basic 
flow was conducted by Hall, [7], He showed that only a very small crossflow velocity is 
needed to radically alter the linear stability properties of short wavelength Gortler vortices. 
In many practical situations interest cannot be legitimately restricted to high wavenumber 
modes alone and, on a linear stability basis, it is frequently the fastest growing (or most 
unstable) vortex which is likely to be important . That this mode is of interest is clear for 
often there will be a whole spectrum of vortex wavenumbers that may be excited in any 
one particular problem and the fastest growing of these is the one that is dominant. The 
most unstable Gortler mode within a two-dimensional boundary layer was identified by 
Denier et. al. [8] and the linear and nonlinear properties of this mode as the crossflow 
component is increased have been developed in [9-11], To the best of our knowledge the 
most unstable vortex in the Taylor problem has yet to be noted. 

The purpose of this concise article is to outline the effect of crossflow on vortices in 
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the Taylor situation for cylinders with a small separation. In particular we shall describe 
the influence of crossflow on high wavenumber modes, specify the structure of the most 
unstable vortex and track the linear stability properties of this mode as the crossflow grows. 
It may be argued that this problem is very specific and therefore is of little, or no, practical 
use. However we should emphasise that the analysis here encompasses the ideas behind 
a large number of important cases including problems involving rapidly rotating bodies 
and impulsive motion. Thus the merit of using the classical small gap Taylor problem as 
a paradigm is that the analysis is very straightforward, enables the important features of 
the flow development to be seen and our workings may be easily generalised to account for 
more complicated situations. 

The remainder of the paper is organised as follows. In §2 we describe the relevant 
stability equations and in §3 show how increasing crossflow modifies the large wavenumber 
structures. The most unstable mode is obtained in §4 and we then outline the influence of 
crossflow on this vortex. We close with a brief discussion. 

2. The governing stability equations 

The linearised disturbance equations pertaining to Taylor vortices have been derived 
many times and so here we restrict ourselves to the briefest of descriptions of these 
equations. Consider the flow of incompressible fluid of density p and kinematic viscosity u 
in the gap between a pair of concentric cylinders of radii R\ and I ? 2 (> R\)- We assume 
that the curvature of the gap is small; that is 8 = (R 2 — R\ )/ R\ "C 1. The inner cylinder 
rotates with angular velocity Uq j R\ whilst the outer cylinder is at rest. It is convenient 
to define variables x,y,z,t by 

Red' r' -Ri z' U 0 t' m 

x_ 8 ’ y ~ R 2 -Ri’ Z Ri-Ri' Re(R 2 -RiY 1 ’ 

where (r',9',z') are the usual cylindrical polar co-ordinates, t 1 is the time variable and 
Re = Uo(R 2 — -Ri)/^ > 1 is the Reynolds number. Writing the velocity field (scaled 
on the characteristic fluid speed U Q ) with respect to ( x,y,z ) as (u,v,w) and scaling the 
pressure p on pUg gives that in the basic state 

(u,v) = (u,v) = (1 - y,0) + 0(6). (2a) 

In order to introduce a crossflow component into the flow we suppose that a constant 
axial pressure gradient acts and the size of this gradient is chosen so that the crossflow 
affects the subsequent stability equations. It is a routine calculation to show that this 
occurs when the axial velocity component w = 0(Re 1 ) and then the remainder of the 
basic state solution is 

(to.p) = (e.p) = A (-^ { 9(1 - y) + 0 («)) , {iz + 0 ( 0 }) ; ( 26 ) 
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where we have introduced A as a measure of the strength of the scaled axial flow. The 
situation discussed in this article has been studied by many prolific authors, the reader is 
referred to [12] for further reading. The earliest theoretical and experimental studies were 
apparently made by [13] and [14] respectively. 

The Taylor number T = 2 Re 2 8 is taken to be an 0(1) parameter and we perturb the 
basic state by writing 


(u,v,w,p) = (u,0,w,p)+A | U(x,y,z,t), j^ v (x, y, z,t), W(x,y,z,t ), -^P(x,y,z,t) j , 

(3) 

where A is a vanishingly small parameter. If we substitute (3) into the linearised forms of 
the continuity and Navier-Stokes equations and neglect terms of relative size 0(6), we are 
left with the four disturbance equations 


dU dV dW_ 
dx dy ^ dz ’ 


( 


d 2 d 2 


d 2 


+ 


dy 2 
d 2 


+ 


dz 2 


d_ 

di 


u=(l - y) W- v+Xy(l - y) ™, 


dy 2 dz 2 


d_ 

dt 


v ~ir = il 

dy 


dV dV 

»> aT + Ai,(1 " y) dl ~ T(1 


( 


d 2 d 2 d\ dP dW 

+ ^-^W-^ = (l-y)^ + \(l-2y)V + \ y (l- y ) 


- y)U , 

dW 

dz ’ 


(4a) 

(46) 

(4c) 


dy 2 ' dz 2 dt, 

which need to be solved subject to the no-slip conditions U = V = W = 0 on y = 0, 1. 
In the absence of crossflow (A = 0) and with no azimuthal dependence (d/dx = 0) system 
(4) reduces to an ordinary differential eigenvalue problem for the neutral Taylor number 
as a function of vortex wavenumber a. It is well known that for aC 1 then T = 0(a~ 2 ) 
whilst, for large a, T = 0(a 4 ). We now analyse the effect of crossflow on this right-hand 
branch of the neutral stability curve; i.e. on neutral large-wavenumber vortex structures. 


3. The effect of crossflow on large wavenumber modes 

It was shown in [1] that large wavenumber vortices in a flow with zero axial component 
are trapped in a thin zone of depth 0(a~3) next to the inner cylinder wall y = 0. 
Analysis of equations (4) reveals that when the scaled crossflow size A = 0(a) this wall 
structure needs to be adjusted and consequently we shall examine this case now. In the 
neighbourhood of the neutral stability curve it is known that the vortices have temporal 
growth rate 0(a 2 ) and a waviness in the azimuthal direction on an 0(a~ % ) lengthscale. 
Using the scalings proposed in [1] we write the Taylor number 

T = T 0 a 4 + ..., 


(5a) 
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seek modes proportional to 


E\ = exp 


laz 


+ a5(o- 0 +a 5 crj 4. . . ,) x 4- a 2 j (fl 0 (f) + a $ (<) + .. -)dt 


( 56 ) 


and allow the crossflow size A = a\. If we further define the 0(1) co-ordinate ( within the 
wall layer by 

y = (6a) 

then the disturbance quantities (U,V,W,P) develop according to 
(U,V,W,P) = 

( Uo + a 5 f/i + . . . , a 2 Vo + a*V\ + ...,a 5 Wo 4~ aW\ a 5 Po 4- (i 2 Pi + . . E\ . 

( 66 ) 

Substitution of (5) & (6) into (46, c) yields at leading orders 


-( l + Cl 0 )Uo = -V 0 , 


-(1 -(- f2o)ko = ~TqUq, 


which are consistent only if 


To — (1 + fi 0 ) 2 - 


(7) 


In order to determine the structure of the leading order disturbance functions it is necessary 
to proceed to higher orders in (4). These give the governing equation 


(ilo + 3)-^^ ho ^1 + fio + 2iX^j — ^—(1 + fto) — 2(f2i -f- <7o) 
which is a scaled version of the Airy equation and has solution 


= 0, 


Vq oc Ai 


1 4“ flo 4" 2* A 
3 4" flo 


t-c 


where 


C = -i Q ° + 3) ! „ ( S-(l + fio) - 2(fi, + <r 0 )) . 


(8a) 


( 86 ) 


(1 4- fio + 2iA)s \To 

In order to satisfy the inviscid constraint that Vo —* 0 as £ — > 0 it is necessary to choose 
C so that Ai(—C) = 0 whence 

C « 2.3381 (8c) 

and then (86) may be rearranged to give the Taylor number correction T \ . Although 
Vo — ► 0 as ( — ► 0 the continuity equation (4a) gives W 0 - i(dV 0 /d£) and so a passive 
viscous layer is required to ensure that the axial disturbance velocity component vanishes 
on the inner cylinder. The details of this layer are very elementary so will not be pursued 
here. 
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Definition (5 6) shows that for neutrally stable large wavenumber modes the frequency 
terms fio', Di are necessarily purely imaginary. Thence result (7) gives To = 1 and (86) 
shows that 

T\ = 3^C(1 + 4A 2 )s cos tan _1 (2A)^ . (9) 

A ~ 2 

Consequently for large scaled crossflow A, T\ = O(As) and expansion (5a) implies that 
when A ~ a the neutral Taylor number is altered by an 0(a 4 ) amount from its value in 
the zero crossflow state. Furthermore, the extent of the vortex activity is compressed into 
a thinner, 0(a -1 ), depth region adjacent to y = 0 and so we write 



X = a 2 X, 

II 

><r 


(10a, 6) 

E 2 = exp 

iaz + a 2 (<7o + 

. . .)x + a 2 J (f2o(0 + • * 


(10c) 

and write the vortex in the form 




{U,V,W,P)=( 

&o + ...,a 2 ^ 0 


• •) E 2 . 

(lOrf) 


Substitution into (4) leads to the disturbance equations 


d 2 

dx 2 



( 11 ) 


which need to be solved subject to U 0 = V 0 = dVo/dx = 0 at x = 0 and as x — ► oo. 

At this point it is worthwhile to recall the results of [ 7 ] which was concerned with 
describing the effect of crossflow on the neutral stability properties of high wavenumber 
Gortler vortices. In the absence of crossflow in that case the mode is confined to a thin 
layer away from the bounding surface. Three distinctive sizes of crossflow were identified. 
In the first the decay of the vortices away from the thin layer is primarily due to the fact 
that the local Gortler number is maximised at the centre of the layer; at a higher crossflow 
the decay is driven by convective effects and in the third regime the layer is driven to 
the wall and the vortex is governed by a pair of viscous equations similar in style to (11) 
but containing quadratic terms as well. Thus in the Gortler (or equivalently the Dean) 
problem the evolution of the vortex with increasing crossflow passes through three phases 
whereas for the Taylor problem studied here we have only two stages. 

Equations (11) were solved using a finite difference technique and, in order to capture 
neutral modes, the quantities Oo and a o were taken as purely imaginary. The results 
expressing To and (f)o + do) as functions of the scaled crossflow A are shown in figure 1. 
Of interest is the solution of (11) in the limits of large or small crossflows. As A — > 0 then 
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we find that T 0 — > 1 and Im + d 0 ) 0 in agreement with the large A results for the 

case presented above. When A — ► oo the solution divides into a double boundary layer 
structure. We postulate that 


Z 2 

Tq ~ dX + . . . , 


flo + do 



+ 


( 12 ) 


for some real constants c and d. In the majority of the region where \ = 0(1) then 


ir 0 = l V.J + ..., V Q = vj + ..., 


where 


(13a) 


-'X (£ 2 - l) V 0 = <®0. 'XPj = V 0 ' ( 13( 0 

and the solution of this system which decays exponentially as y — > oo is 

WW) = (x-^KvW.ix^KM) , (M) 

where K u is the modified Bessel function of order 1 / = yjd+ j (see [15]). As x — * 0 these 

^ — I 

ev. 3 ~ 

solutions develop singularities so that in a wall layer where, x = A Y,Y = 0(1), we write 

(15a, b) 


ir 0 = P °u}(Y) + 




These functions satisfy the coupled equations 


( 


d 2 

dY 2 



•pyj 

dY 2 





u} = -v; 


0 i 


(15c, d) 


+ + + ~ ~ 

which need to be solved subject to the conditions that U 0 = Vq — dV 0 j dY = 0 on Y = 0 
together with suitable matching conditions to (14) as Y — > oo. Numerical solution of this 
problem reveals that 

c = -4.706, d = 0.721, (16) 

and the asymptotic forms (12) for To and + do are sketched on figure 1. As may be 
seen, agreement between the numerical solution of (11) and these asymptotes is eminently 
satisfactory and demonstrates that as the (scaled) crossflow moves through this 0(a 2 ) 
regime the changes induced in the neutral Taylor number are substantial. Therefore only 
a very small axial field of size 0(i?e _1 a 2 ) is required in order to radically distort the 
neutral curve for large wavenumber modes and so suggests that any such wavenumber 
vortex would be completely stabilised by a tiny crossflow. 

In many practical situations Taylor/ Gortler numbers are large and, whilst the effect of 
any crossflow component on the neutral curve is of interest, it is also important to describe 
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how the crossflow modifies the growth rate of the most unstable vortex. As was mentioned 
in the introduction, the most unstable mode is frequently the one that is dominant in 
practice and, as far as we are aware, the location of this mode has not been written down 
for the Taylor problem. Therefore in the coming section we first obtain the most unstable 
vortex in the absence of crossflow and then consider how the imposition of a gradually 
increasing crossflow component modifies the stability properties of this disturbance. 

4. The most unstable mode 

In order to deduce the most unstable linear vortex in the absence of crossflow it is 
easiest to proceed in the manner of Denier et. al. [8]. In this work the authors considered 
both the growth rate of-inviscid (0(1) wavenumber) Gortler modes as the wavenumber 
increases and also the properties of modes in the vicinity of the right hand branch of the 
linear neutral curve. It was then possible to identify an intermediate wavenumber regime 
in which infinitesimal high Gortler number vortices are the fastest growing. 

The inviscid problem relevant to Taylor vortices for T 1 has been considered 
numerous times (see [16] for example) and for vortex wavenumber a = 0(1) the 
disturbance takes the form 
(U,V,W,P) = 

(U 0 (y) + . . . , TW 0 (y) + . . . , T?W 0 (y) + . . . , TP 0 (y ) + . . .) exp jiaz + T^ J «(*)<**} • 

Substitution of these forms into the governing equations (4) shows that Vo satisfies 

ft 2 - a 2 V 0 ) = -a 2 (l - y)V 0 , (17) 

subject to the boundary conditions that Vo = 0 at y = 0, 1. This gives an eigenproblem 
for the growth rate as a function of a. In order to locate the most unstable vortex across 
the entire wavenumber spectrum it is useful to consider the behaviour of the solution to 
(17) as a — > oo. Elementary analysis (see [2]) shows that the eigenfunction Vo is confined 
to a thin zone of depth 0(a~ '^ ) attached to y = 0. Thence 

n = 1 - \CaT$ +..., (18) 

where the constant C satisfies Ai(-C) = 0 and has been given in (8c). Thus as a — * oo 
the growth rate of the vortex asymptotes 

T^l-lCa-i + ...). (19) 

Otto [2] obtained this result for his problem concerning the impulsive spin-up of a cylinder 
within an expanse of fluid and we conclude that for large a the growth rate becomes 
independent of wavenumber to leading order. 
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Results (5a) and (7) demonstrate that at large vortex wavenumbers the right hand 
branch of the neutral stability curve is given by T = a 4 + — Then for T 1 if we 
examine a vortex of wavenumber a — aT « + . . . with a = 0(1) it is straightforward to 
deduce that the vortex is confined to a region of depth 0(T~ £ ) next to the inner cylinder 
wall. Furthermore, the sizes of the disturbance components (U,V,W, P) are in the ratios 
(1, T^, , Ts) and the leading order vortex growth rate term is 0(T%). If we let this 
growth rate Cl = T^Cl + . . . then the compatibility criterion arising from the leading order 
balances from momentum equations (46, c) reduces to 

Cl = 1 — a 2 . 


Thus, as would be expected, the vortices are unstable for a < 1 and stable for a > 1. In 
addition, the overall growth rate is 

T? (l — a 2 ) + . . . (20) 

Mid thus the expressions (19) and (20) become comparable when 



a ~ 


TTe, 


3 

This suggests that the most unstable Taylor vortex has wavenumber 0(T T*) and 
we shall confirm this presently. Meanwhile we note how the geometry of growth rate/ 
wavenumber space for the Taylor problem contrasts with that for the Gortler (or Dean) 
case. We have seen from (19), (20) that for a large wavenumber range, specifically 
1 < a < T?, the growth rate of the Taylor mode is constant at leading order and only 
varies in higher order terms. Conversely, in the Gortler case the leading order growth rate 
increases both as a — » oo in the inviscid problem and as a — > 0 in the right-hand branch 
case. Therefore in the intermediate regime where the growth rate is greatest the maximum 
is a much more pronounced feature. This difference in growth rate behaviours for the two 
problems does have important consequences for the subsequent analysis. For the Gortler 
case the calculations of Bassom & Hall [9] show that steadily increasing crossflow stabilises 
the vortices and only a small crossflow is required in order to ensure that the flow is stable 
to infinitesimal modes irrespective of their wavelength. In addition it was found that 
as the crossflow grows so the location of the most unstable vortex remains in the same 
wavenumber regime: behaviour which is not repeated for the Taylor case discussed below. 

In order to prove that it is the 0(T A) wavenumber regime which contains the most 
unstable Taylor mode we write 

a = a 0 T^, (21a) 


and seek disturbances of the form 


(U,V,W,P)= (Uo + ...,T±Vo + 


.,T*Wo + ...,TfPo + ...)^3, (216) 
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where 


E 3 = exp 


ia 0 T&z + T?J (ft 0 (*) + T"* 6, (*) + -..) <ft] , (21c) 


and Uo,Vo, . . ■ are functions of i) = T»y. Therefore these vortices axe trapped in an 0(T « ) 
thick zone adjacent to y — 0 and we shall impose a crossflow of scaled size 


X = Tre\. 


( 22 ) 


This scaling is arrived at by choosing sufficient crossflow so as to perturb the vortex 
structure relevant to the zero crossflow problem without grossly changing its underlying 
characteristics. We substitute (21), (22) into equations (4) and compare like powers of T. 
The consistency condition arising from momentum balances (46, c) leads to 




= 1 , 


(23) 


as would be anticipated from result (20). At next order we find that Vo satisfies the scaled 
Airy equation 



di 7 2 


— Gq(1 T 2i\ao) 



2(ajj T fli ) 

(IT 2iAdo) 



Demanding that Vo = 0 at r] = 0 and as r] —> oo leads to an expression for Di and thence 
we deduce the overall vortex growth rate Gr = T 5 T £l\T% T o(Ts) where 


D] = 


C{ 1 T 4doA 2 )3 


2a 0 3 


|c.os tan ! (2doA)^ + i sin tan 1 (2«oA)^| — dg. (24) 


A - _ £ * 

In the absence of the crossflow component Re{Q, i) = — ^Cd 0 3 — Oq which clearly —* — oo 
as do — ► 0 or as do — * oo. The maximum value of this expression is — 4((7/6)* ~ —1.973 
which occurs when do = (( 7 / 6)5 « 0.702 and consequently we have the result that in the 

3 

classical Taylor problem the the most unstable mode has wavenumber ~ 0.702 Tt«. As the 

A A | A A 2 

scaled crossflow increases then Re( fli) is greatest when do = 0( A“») and |fli| = 0(A3). 
Furthermore, for a fixed do, as A increases so the overall vortex growth rate diminishes and 
crossflow is seen to have a stabilising effect. It is clear that a revised solution structure 
is needed when the second term in Gr becomes comparable with the first, i.e. when 
= 0(T?) or A = 0(Trs). 

We therefore define 

A = T*A, (25a) 

and seek the wavenumber of the most unstable vortex for this case. We find that for a large 
range of wavenumbers, specifically l<a<C , the leading order component of the normal 

velocity field, say Vo-, is given by the solution of a Whittaker-type equation (see [15]). The 
disturbance structure is confined to lie within a thin zone adjacent to y = 0 whose precise 
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depth is a function of the wavenumber. However for all vortex wavenumbers 

the leading order disturbance forms are determined by the solution of a common equation 

which may be solved analytically in terms of the modified Bessel function K v . In order to 

satisfy both the no-slip conditions at the wall and suitable decay conditions at the edge of 

* 1 

the thin zone it is found that the leading order growth rate of the vortices is T* where 
fio is a solution of the equation 

Kp(ho/i\)=0 t ji 2 = \ + i 2 . (25 6) 

ys - - £ 

The dependence of CIq upon A is shown in Figure 2. When A = 0 vortex modes with 
wavenumber 1 ^ a T* have growth rate o = 1. As A increases so the growth 
rate diminishes. Indeed, as A — * l/\/2 so — * 0 — i/\/2 and all the modes with 

wavenumber within 1 a T 4 are made neutral. Again therefore, just as in the Gortler 
problem, crossflow is seen to have a stabilising influence and at crossflows of greater than 
Re~ l T % /\/2 the flow configuration is stable to vortices of any wavenumber CTc 

One question we have not addressed here is that concerning the structure of these vortex 
modes as the crossflow parameter A increases through the critical value l/\/2- Indeed 
when A = I/a/ 2 precisely then the formal solution for F 0 develops a singularity and the 
disturbance structure requires a critical layer type zone in order to smooth this out. This 
situation has similarities with the scenario discussed by Blackaby & Choudhari [17] who 
were concerned with inviscid vortex modes in slightly three-dimensional boundary layers 
and how such perturbations are related to the Taylor-Goldstein equation for stratified shear 
flows. These authors found solutions for which the leading order ‘vertical’ (i.e. normal to 
wall) velocity component of the inviscid disturbance becomes infinite at the critical layer 
and they asserted that theirs was the first case in which this type of behaviour has been 
found. We also believe this to be the case and suggest that we have here another problem 
where this phenomenon is seen. 

The discussion above implies that for a large range of vortex wavenumbers all modes 
have the same leading order growth rate. However, such vortices axe not the fastest growing 
ones as further investigation reveals that within the 0(T?) wavenumber regime we have 
the unique most unstable mode. In order to demonstrate this we define the crossflow 
parameter A as in (25a) and look for disturbances of wavenumber a = aT? . If we let 

F /4 = exp j^aT 4 2 T* J ^fi(t) + . . dt ^ , (26a) 


and seek solutions of the type 

(UyVyWyP) = ( U + . . . , T ^ V + . . • , T ^ W + . . . , T ^ P + . . . )E 4 , 


( 266 ) 
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where U, V, W, P are functions of the scaled co-ordinate </> = T * y we find that U and V 
axe given by the coupled viscous equations 


d? 
d<j > 2 



- 1 JV = a 4 U , 

u = -v, 


(27 a) 


(276) 


subject to the conditions that U = V = dV /d<f> = 0 at <j) = 0 and as (f> —*■ oo. As would be 
expected in view of the scalings, these equations are essentially identical to those relevant 
to determining the right-hand branch of the neutral curve, (11). 

We solved (27) for the (in general complex- valued) growth term as a function of the 
crossflow parameter A and scaled wavenumber a, and the results of these calculations are 
summarised in figure 3. This figure shows that for A < 1 / y/2 the vortex is unstable over 
a finite range of wavenumbers including those with d — > 0 and the extent of this unstable 
range decreases with increasing A. Indeed when A reaches the value 1 / \/2 the mode with 
d = 0 is made neutral in agreement with the results of the analysis relevant to modes 
with wavenumber 1 <«< . For increasing A beyond l/\/2 the unstable band reduces 

further until when A ~ 0.942 the whole wavenumber spectrum is made stable. 

The results given in figure 3 also show the dependence of the most unstable vortex 
wavenumber upon the crossflow A. If we denote the wavenumber of this most unstable 
mode by d u then as A — > 0 so d u — > 0 since we know from our previous workings that as 
the crossflow diminishes so the most unstable vortex reverts to within the 0(Ti«) regime. 
In addition, d u is not a monotonic function of A and as A — 0.942 it is the mode with 
wavenumber d ~ 0.103 which is the last to be made stable. 

The calculations in this and the preceding section give related but not identical 
information. The work of §3 may be regarded as describing the stability properties of 
a vortex of specified, large wavenumber. In particular, the solution of equations (11) 
accounts for the dependence of the neutral Taylor number on the size of the imposed 
crossflow whereas the work of this section describes how crossflow affects the stability of 
vortices within a flow at a specified large Taylor number. We notice how the most unstable 
vortex has wavenumber O(T^) in the absence of crossflow; with increasing crossflow this 
fastest growing mode moves into the O(T^) wavenumber regime. At this stage the right 
hand neutral branch of the stability curve is also within the O(T^) range so that we can see 
that as the crossflow component increases from zero so the location of the most unstable 
linearised mode moves out into the asymptotic regime containing the right hand branch. 
This contrasts with the. situation in the Gortler case where, for a large Gortler number G 
and in the absence of crossflow the right hand branch is at an 0(G^ ) wavenumber whilst 
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the fastest growing mode has wavenumber 0(0* ), [8], As the crossflow increases the 
most unstable vortex remains within the 0 ((j 5 ) regime and the right hand neutral branch 
moves into this regime so that it is a mode of wavenumber O(G^) which is the last to be 
stabilised. 


5. Conclusions 

In this work we have studied how the imposition of a small crossflow velocity 
component affects the linear stability characteristics of vortices in high Taylor number 
flow. In particular we have identified the most unstable vortex, illustrated ihat for 
large wavenumbers a crossflow of size 0(Re is sufficient to completely stabilise 

the mode and proved that crossflow of 0(Re~ 1 T ?) makes the flow immune to unstable 
vortices irrespective of their wavenumber. In some ways it is this last result which is 
the more significant for it shows that in many flows where an appreciable degree of 
three -dimensionality is present the Taylor vortex mechanism is likely to be unimportant 
compared with other instability forms. 

Our results may be compared with those of the equivalent Gortler or Dean problem 
in the following way. Suppose that instead of having the inner cylinder of our geometry 
moving with speed Uq it is at rest and an azimuthal pressure gradient is imposed of 
sufficient strength so as to induce an 0{U Q ) fluid speed across 0 < y < 1. The analysis 
of [8] suggests that in the absence of crossflow the most unstable vortex has wavenumber 
0(T*). In addition a crossflow of size 0(J?e -1 T t ) would be needed to completely stabilise 
the flow— a value which is seen to be greater than that in the Taylor problem considered 
here. This suggests that in situations which are composed partly of a Taylor-type situation 
and partly a Dean-type then it the latter component which is dominant in determining the 
stability characteristics when crossflow is present. However, there will still be circumstances 
in which the structure presented in this paper is important and such cases include spin-up 
configurations of the type studied by Otto [2]. In that problem a cylinder at rest in an 
incompressible fluid is given an impulsive angular velocity. A Rayleigh layer is set up 
and this flow is prone to vortices which are identical in character to those examined here. 
If a basic flow component is added along the axis of the cylinder (so this could model 
a spinning cylindrical object moving through a fluid) this axial flow may be determined 
using the method of Glauert & Lighthill, [18]. 

This applicability of our analysis to other flow problems is an attraction of our working. 
The flow configuration chosen at the outset was deliberately taken to be very simple but 
it does capture most of the essential features of a wide class of flows over rotating bodies 
whenever the inviscid Rayleigh criterion for instability is most violated at the body surface. 
Not only can our analysis be extended to other geometries but it may also be easily adapted 
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to allow for vortex motions modulated in the azimuthal direction. Consequently it has the 
potential of proving to be useful in a number of related problems. 

There are several directions in which our calculations might be extended. The most 
obvious of these is the question of nonlinearity. In order to give a definitive account of 
crossflow effects it would be necessary to allow for other than infinitesimal vortex modes 
by pursuing weakly nonlinear or strongly nonlinear computations. Equivalent work for 
the Gortler problem is reported in [10, 11, 19] where it is shown that weak nonlinearity 
stabilises near neutral but linearly growing modes to give non-zero finite amplitude vortex 
states. Inclusion of full nonlinearity tends to suggest that the vortices break down in a 
finite distance singularity which would be the analogue of a finite time singularity in the 
present work. In addition, of theoretical interest would be a rigorous analysis of the critical 
layer structure which arises just as vortices of wavenumber l<a<cT« near neutrality as 
the crossflow approaches Re~ l Ti j\J 2. This analysis would certainly be required in order 
to investigate nonlinear configurations near this cut-off. However the practical significance 
of such considerations is far from clear at this juncture for at this crossflow size we have 
shown that there is still a wide range of higher wavenumber vortices which are linearly 
unstable. Indeed it is not until that the axial flow becomes greater than ~ 0.942Z?,e~ 1 
that infinitesimal Taylor vortices of arbitrary wavenumber are made stable. 

We are grateful to Dr. J.P. Denier for his valuable comments on an earlier draft of this 
work. The research of the first author was supported by the National Aeronautics and 
Space Administration under NASA contract No. NASl-19480 while he was in residence 
at the Institute for Computer Applications in Science and Engineering (ICASE), NASA 
Langley Research Center, Hampton, VA 23681, USA. 
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Figure 1 . Solution of equations (11) in order to determine the location of the right-hand 
neutral branch as a function of crossflow parameter A; a) T 0 ; b) Im(f2 0 + <r 0 ). The dotted 
lines on the figures denote the asymptotic forms (12) and (16). 
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Figure 3. Solution of viscous equations (27) to determine leading order vortex growth 
rate Re(Q)T? as a function of vortex wavenumber aT* and crossflow A . The graphs 
correspond to crossflow values A = 0.02, 0.05, 0.1 mid thereafter at intervals of 0.1 up to 
and including A = 1.1. The A symbols indicate the locus of the most unstable vortex 
wavenumber. 
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